ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY 
COUPLED SYSTEMS OF PDE 
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Abstract. The adjoint method, recently introduced by Evans, is used to study obstacle prob- 
lems, weakly coupled systems, cell problems for weakly coupled systems of Hamilton-Jacobi 
equations, and weakly coupled systems of obstacle type. In particular, new results about the 
f ^ , speed of convergence of common approximation procedures are derived. 
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1. Introduction 

In this paper we study the speed of convergence of certain approximations for obstacle problems 



Pl' 1 an d weakly coupled systems of Hamilton-Jacobi equations, using the Adjoint Method. This 

technique, recently introduced by Evans (see jEvalOj . and also |Tra] and |CGTj ). is a very successful 



tool to understand several types of degenerate PDEs. It can be applied, for instance, to Hamilton- 
Jacobi equations with non convex Hamiltonians, e.g. time dependent (see [Eva 10] ) and time 
independent (see |Tra ) Hamilton-Jacobi equations, to weak KAM theory (see |Fat97ai [Fat97bi 
IFat98a[|Fat98b] . [EGOlj . |EG02] ). and to the infinity Laplacian equation (see |ESj ). We address 
here several applications and propose some new open questions. Further results, which will not 



> 

OA , 

■ be discussed here, can be found in [Eva 10] and |CGT| . 

cn 

1.1. Overview of the Adjoint Method. To apply the Adjoint Method to a (non linear) PDE, 
one has to consider the adjoint equation associated to the linearization of the original problem. 
In this way it is possible to prove new estimates, which can then be used to obtain additional 
information on the solution of the initial PDE. In order to give an idea of the technique, we 
show below how this was used in the context of Aubry-Mather theory, in the periodic setting (see 
[CGT] V 

To start with, we quote a fundamental result (see [LPV88] ). stating existence and uniqueness 
of the effective Hamiltonian. Here with T" we denote the n-dimensional unit torus in K™, n £ N. 



Theorem 1.1 (Lions, Papanicolaou and Varadhan) . Let H : Fxl" — > K be smooth and coercive, 
i.e. 

lim H(x,p) = +oo. 

|p|— >+oo 
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Then, for every P £ R™ there exists a unique H[P) £ R such that the equation 

H(x, P + D x u(x, P)) = H(P) (1.1) 

admits a Z n -periodic viscosity solution u(-,P) : T n — > R. 

In the spirit of Theorem ll.il we prove in [CGT] the analogue result for an elliptic regularization 
of equation (JTTTJ). See also |Gom02| for similar results. 

Theorem 1.2 (See [CGT] ). Let H : T" x R™ -> R be smooth and assume that 
lim (-\H(x,p)\ 2 + D x H(x,p) -p) = lim H ^f^ = +00 . 



|p|— >-+oo 



2 J \pH+oo \p\ 



Then, for every r/ > and euen/ P <G R™, t/iere exists a unique number H* 1 (P) £ R sitc/i t/iat </ie 
egitation 

- ^ Au"(i) + JT(z, P + Du v (x)) = W{P) (1.2) 
admits a unique (up to constants) Z n -periodic viscosity solution. Moreover, for every P £ R" 

lim H V (P) = H(P), (up to subsequences) 

where H(P) is given by Theorem ] 1. 11 In addition, we have 

u n — > u uniformly, (up to subsequences) 

where u is a Z" -periodic viscosity solution of (|1.1[) . 

For every -q > and P £ R ra , the formal linearized operator L r,,p : C 2 (T") — > C(T n ) associated 
to equation (|1.2p is defined as 

2 

L"< p w (x) := _1-A«(ar) + D p H(x, P + DvP(x)) ■ Dv{x), v £ C 2 (TP). 

As already mentioned, the main idea of the method consists in the introduction of the adjoint 
equation associated to L n,p : 

2 

-?-Aa n (x)-dw{D p H{x,P + Du v {x))a v (x)) = 0, in TP, 

, (1-3) 
a 11 is non- negative, T"-periodic and / d n dx = \. 

Then, exploiting the properties of the solution a 71 of (|1.3[) . we can retrieve additional information 
about u v . 
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First of all, one can show new estimates that do not seem to be easily obtained in a classical 

\Du r '\ 2 

way. As an example, define the function «/' := — - — . Then, w 11 satisfies 

2 2 

D X H ■ Du'i + D p H ■ Dw 1 ! = yAw" - y |£>V| 2 . (1.4) 
Multiplying the above relation by and integrating by parts, we eventually get 



if 



|£»V|Vdx< C, (1.5) 



for some C > independent of r\. Relation (|1.5[) gives information about the behavior of all 
the Hessian D 2 u v of u v in the support of a n . We observe that, without passing to the adjoint 
equation, one can only conclude that 

7? 2 |Am"| < C, 

thus obtaining a relation which involves just the Laplacian AvP of u v . More generally, by consider- 
ing functions of the form w v (x) = 4>{x, P + Du n (x)) and studying the analogous of equation fj 1 .4f) . 
one can obtain further properties, using compensated compactness based estimates (see [Eva 10] . 

[CUT]). 

In addition, the Hamiltonian H is not required to be convex. When we have that H is uniformly 
convex in p (i.e. D 2 p H > a for some a > 0), (|1.5p can be significantly improved (see [CGTj V 
Indeed, differentiating (|1.2I) twice along a generic vector £ £ M. n with |£| = 1 we have (here and 
always in the sequel, we use Einstein's convention for repeated indices in a sum) 

2 

ff« + 2fl c *4, 1 + H PiPj ul Xi ul Xj + D P H ■ Dv?x = yAu| ? . 
Then, multiplying by a v and integrating by parts we get 



/ \Du'l\ 2 a" dx <C, 



or more generally 

l^i/ftr' dx < C, 



which is clearly stronger than fj 1 . 5f) . The differences between convex and nonconvex setting can 
be also observed by investigating the existence of invariant Mather measures (see |CGT| V and by 
studying the nature of the shocks in Hamilton- J acobi equations (see [EvalOp . 

Finally, the treatment allows to analyze the speed of convergence oiH'(P) to H{P). Indeed, 
classical arguments in elliptic regularization imply that u n and H v are smooth in r\ away from 
rj = 0. Then, differentiating (jl.2l) w.r.t. r\ 

DpH ■ Dv% = H^iP) + \&vP n + V Au\ 
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where we denoted the differentiation w.r.t. rj with a subscript. Again, multiplying by a 71 and 
integrating by parts we infer that 

\H^(P)\<riJ \Au ri \a' n dx<Cr](J |I>V|V»cteV (j o* dx\ ' < C, 
where the latter inequality follows by using (jl.5p . Thus, we conclude that 

\TT{P)-H(P)\<c^ 

which shows that the speed of convergence is 0(rj). In the uniformly convex case, this result can be 
improved to 0(r] 2 ). Let us observe that, as pointed out in [Eva 10] and by Fraydoun Rezakhanlou 
to us, these estimates on the speed of convergence can also be obtained via Maximum Principle. 

1.2. Outline of the paper. This paper contains four further sections concerning obstacle prob- 
lems, weakly coupled systems, effective Hamiltonian for weakly coupled systems of Hamilton- 
Jacobi equations, and weakly coupled systems of obstacle type, respectively. We use a common 
strategy to study all these problems. For this reason, we describe in more detail our approach just 
in the case of the obstacle problem. In all the paper, U is an open bounded domain in R" with 
smooth boundary, n > 2. Moreover, we will denote with v the outer unit normal to dU. 

In Section [2] we consider an obstacle problem of the form 

max{« — tp, u + H(x, Du)} = in U, 

(1.6) 

u = on dU, 

where ip : U — >• R and H : W l x U — >• R are smooth, with -0 > on dU. This equation arises 
naturally in Optimal Control theory, in the study of optimal stopping (see |Lio82j l In |IY90] . 
Ishii and Yamada consider similar problems in the special case in which H is uniformly convex in 
the second variable. 

Classically, in order to study (I1.6P one first modifies the equation, by adding a perturbation 
term that penalizes the region where u > tp. Then, a solution is obtained as a limit of the solutions 
of the penalized problems. More precisely, let 7 : R — > [0, +00) be smooth, such that 

7(s) = for s < 0, 7(5) > for s > 0, 
< j'(s) < 1 for s > 0, and lim j(s) = +00, 

and define j £ : R — > [0, +00) as 

7 e (s) := 7 ^-^ , for all s € R, for all e > 0. (1.7) 
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In some of the problems we discuss we also require 7 to be convex in order to obtain improved 
results, but that will be pointed out where necessary. For every e > 0, one can introduce the 
penalized PDE 

u £ + H(x, Du £ ) + j £ (u £ -ip) = eAu £ in U, 

(1.8) 

u £ = on dU. 

To avoid confusion, we stress the fact that here 7 e (u 6 — ip) stands for the composition of the 
function j £ with u £ ' — tp. Unless otherwise stated, we will often simply write 7 e and (7 e )' to denote 
7 e (u e — ip) and (7 e )' (u £ — ip), respectively. Also, notice that in (ll.8[) the parameter e corresponds 
to rj 2 in (|1.2[) . We made this choice in order to compare our results with existing estimates for the 
speed of convergence in literature. 

Thanks to [Eio82 , for every e > there exists a smooth solution u £ to (|1.8p . It is also well 
known that, up to subsequences, u £ converges uniformly to a viscosity solution u of Q1.6P (see also 
Section [2] for further details). 

In [IY90j Ishii and Yamada considered related problems when H (x, ■) is uniformly convex, and 
studied the speed of convergence of the functions u £ to u. However, both the original problem, 
the regularized PDE, and their methods are different from ours. To the best of our knowledge, no 
results are available in literature concerning non convex Hamiltonians. 

We face here the problem requiring a coercivity assumption on H and a compatibility condition 
for equation (|1.6[) (see hypotheses (HUl) an d (H2J2), respectively), showing that the speed of 
convergence in the general case is 0(e 1 ^ 2 ). 

Theorem 1.3. Suppose conditions an d (1^2) in Section^ hold. Then, there exists a 

positive constant C , independent of e, such that 



U - U\\ L a 



< Ce 1 ' 2 . 



At the end of the section, we give a dynamic and a stochastic interpretation of the problem 
(see Subsection 12.31 and Subsection 12. 4[ respectively). The proof of Theorem 11.31 consists of three 
steps. 

Step I: Preliminary Estimates. We first show that 

u £ (x) - V>(s) 

max < C, (1.9) 

for some constant C > independent of e (see Lemma r2"?2"j) . This allows us to prove that 

\\u £ \\ L ~,\\Du £ \\ L ~ <C, 

see Proposition 12. II 
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Step II: Adjoint Method. We consider the formal linearization of fj 1 .8|) . and then introduce 
the correspondent adjoint equation (see equation (|2.6|1 ). The study of this last equation for different 
values of the data allows us to obtain several useful estimates (see Lemma |2~31 and Lemma [53]). 

Step III: Conclusion. We conclude the proof of Theorem 11.31 by showing that 

C du e 
maxK(z)|<™ u%{x) ;= — (as), (1.10) 

xeu £ 1 os 

for some constant C > independent of e (see Lemma T2.5p . The most delicate part of the proof 
of (jl.lpp consists in controlling the term (see relation (|2.11[l ) 

7|(*)=-^7'Q, s€R, 7 |( S ) : =^( S ). (1.11) 



We underline that getting a bound for (|1.11[) can be extremely hard in general. In this context, 
this is achieved by differentiating equation (|1.8p w.r.t. e (see equation (|2.10p ). and then by using 
relation (jl.9l) . Lemma l2.3l and Lemma f2.4l This means that we overcome the problem by essentially 
using the Maximum Principle and the monotonicity of 7 6 (see estimates (|2.12p and (|2.13p 'l . We 
were not able to obtain such a bound when dealing with homogenization or singular perturbation, 
where also similar terms appear. We believe it would be very interesting to find the correct way 
to apply the Adjoint Method in these situations. 

In Section [3] we study the weakly coupled system of Hamilton- Jacobi equations 

cnMi + C12U2 + Hi(x, Dm) = 0, 

in U. (1.12) 

C21M1 + C22U2 + H 2 (x, Du 2 ) = 0, 

Under some coupling assumptions on the constants (see conditions (102) and (103)), Engler and 
Lenhart [EL91j . Ishii and Koike [IK9T] prove existence, uniqueness and stability for the viscosity 
solutions (ui, U2) of (|1.12p . but they do not consider any approximation of the system. 

As before, we introduce perturbed problems (see (13.21) ) and show that, under the same as- 
sumptions of [E L91j . the speed of convergence of the corresponding solutions (uf,u|) to (iti,M2) 
is 0(e 1 / 2 ) (see Theorem 13. 6p . We observe that the coupling assumptions here play a crucial role 
and cannot be replaced. For the sake of simplicity, we just focus on a system of two equations, 
but the general case can be treated in a similar way. 

Section H] is devoted to an analog of the effective Hamiltonian problem (jl.ip introduced by 
Lions, Papanicolaou and Varadhan [LPV88 , which is the following weakly coupled system of 
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Hamilton- Jacobi equations: 

c\U\ — c\U2 + Hi(x, Du\) = Hi 

in T™. (1.13) 

-c 2 ui + c 2 u 2 + H 2 (x,Dui) = H 2 
Here ci and c 2 are positive constants and Hi,H 2 : T™ x R™ — > R are smooth, while Ux, u 2 : T n — > R 
and H\,H 2 £ R are unknowns. Systems of this type have been studied by Camilli, Loreti and 
Yamada in 1CL08) and |CLY09] . for uniformly convex Hamiltonians in a bounded domain. They 
arise naturally in optimal control and in large deviation theory for random evolution processes. 

As in |CL08j . we assume the system to be quasi-monotone and not necessarily monotone (see 
condition (1212)). In this context, by monotonicity we mean exactly the coupling assumptions 
of Engler and Lenhart in |EL91j . Moreover, we also require Hi, H 2 to be coercive (see condition 
(Hill))' As it happens for the cell problem in the framework of weak KAM theory, there is no 
hope here of a general uniqueness result for (ui, u 2 ), even modulo the addition of constants. 

In the spirit of Theorem ll.il studying a perturbation of (I1.13[) (see (I4.2[) ) we prove that there 
exist Hi, H 2 £ R such that the system above admits viscosity solutions ui,u 2 . Notice that Hi, H 2 
are not unique as well, but c 2 Hi + ciH 2 is. More precisely, we show that there exists a unique 
\x £ R such that 

c 2 Hi + ciH 2 = ll, 

for every pair (Hi,H 2 ) £ R 2 for which viscosity solutions ui,u 2 exist (see Theorem 14. 2\ . This 
result is the analog of Theorem II .11 by Lions, Papanicolaou and Varadhan in the case of systems, 
Finally, we prove that the speed of convergence of (eu\,eu 2 ) to (Hi,H 2 ) is 0(e) (see Theorem 
14. 4[) . where itf and u 2 are the solutions of the approximating problem (|4.2p . 

In Section [5j we conclude the paper with the study of weakly coupled systems of obstacle type, 
namely 

max{ui — u 2 — ipi, ui 4- H\(x, Dui)} — in U, 

(1.14) 

max{w2 — U\ — ip2 7 v>2 + #2(^7 -£^2)} = in U. 
Problems of this type appeared in |CDE84] and |CLY09j . Here H x , H 2 : U X R" -> R and 
ipi, "4>2 '■ U — > R are smooth, with ipi, ip 2 > a > 0. 

In this case, although the two equations in (11 . 14)) are coupled just through the difference ui — u 2 
(weakly coupled system), the problem turns out to be considerably more difficult than the corre- 
sponding scalar equation f| 1 . 6[) . Indeed, we cannot show now the analogous of estimate fj 1 .9[) as in 
Section H 

For this reason, the hypotheses we require are stronger than in the scalar case. Together with 
the usual hypotheses of coercivity and compatibility (see conditions (Ej5j2) and (H5]4)), we have 
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to assume that Hi(x,-) and Hzix,-) are convex (see (E(5jl)), and we also ask that D X H\ and 
D X H2 are bounded (see (E(5J3)). We were not able to relax these conditions. We believe it would 
be also very interesting to apply the Adjoint Method in this particular system with the above 
conditions relaxed. 

2. Obstacle problem 

In this section, we study the following obstacle problem 

max{u — ip, u + H (x, Du)} = in U 

(2.1) 

u = on dU, 

where ip : U R and H : W l x U R are smooth, with ip > on dU. We also assume that 

(lEll) lim [ -\H{x,p)\ 2 + D x H(x,p) ■ p] = lim H( ff^ = +00 uniformly in x eU; 

H-H-oo \2 ' J |pK+oo \p\ 

(IEJ2) there exists a function $ e C 2 ([7) n C X (JJ) such that $ < ip on [7, $ = on dU and 

$ + H{x, L>$) < in U. 



We observe that in the classical case H(x,p) = H(p) + V(x) with 

hm — — — = +oo, 

or when iJ is superlinear in p and |D x ff < C(l + |p|), then we immediately have (tJ2]l). 
Assumption (EJU2) (stating, in particular, that $ is a sub-solution of ^.lp ). will be used to derive 
the existence of solutions of (12. ip , and to give a uniform bound for the gradient of solutions of the 
penalized equation below. 

2.1. The classical approach. For every e > 0, the penalized PDE is the equation given by 

u £ + H(x,Du £ ) +j £ (u £ - ip) = eAu e in U, 

(2.2) 

u £ = on dU, 

where 7 e is defined by (11.71) . From [Lio82 it follows that under conditions (E|2]l) an d (E02), for 
every e > there exists a smooth solution u £ to (|2.2[) . The first result we establish is a uniform 
bound for the C^-norm of the sequence {u £ }. 

Proposition 2.1. There exists a positive constant C, independent of e, such that 

\\u £ \\ L ~,\\Du £ \\ L ~ < C. 

In order to prove Proposition ^. 11 we need the following fundamental lemma: 
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Lemma 2.2. There exists a constant C > 0, independent of e, such that 

u £ — ib 

max7 e (u e — ip) < C, max < C. 

xeu xeu £ 

Proof. We only need to show that max xG jj~f £ (u £ — ip)<C, since then the second estimate follows 
directly by the definition of 7 e . Since u £ — ip < on dU , we have max^ga^ 7 e (m £ — tp) = 0. 

Now, if max x( -jjj E (u E — ip) — 0, then we are done. Thus, let us assume that there exists 
x\ G U such that max je jj 7 e (u £ — tp) = 7 e (u e — ip)(xx) > 0. Since 7 e is increasing, we also have 
max x g(/(u E — ip) = u E (xi) — ip(x\). Thus, using (|2.2[) . by the Maximum principle 

(u £ (xi) - ip{x x )) + i £ (u £ ( Xl ) - ip{ Xl )) = eAu £ (xi) - H{x u Du £ { Xl )) - ip( Xl ) 
< eA^(xi) - H(xi, Dip[ Xl )) - ip{x x ). 

Since u £ {x\) — ip{x\) > 0, 

7 e (u e (x 1 ) - xP( Xl )) < max(|AV| + \H(x,D^)\ + \iP(x)\) < C, 
xeu 

for any e < 1, and this concludes the proof. □ 

Proof of Provosition \2. 1\ Suppose there exists Xq € U such that u £ (xq) = max^g-jj u £ (x). Then, 
since Au £ (xo) < and using the fact that j £ > 

u £ (x ) = eAu £ (x ) - H(x ,0) - j £ (u £ (x ) - ip(x )) 

< -H(x Q ,0) < max (-iT(ai, 0)) < C. 

xeu 

Let now x\ € U be such that u £ (x\) = min x£ -jju £ (xi). Then, using Lemma 12.21 
w e (iri) = e&u £ {xi) - H(x u 0) - 7 £ (u £ (xi) - ip(xx)) 

> -H(x 1 ,Q)-'f(u e (x 1 )-il>(x 1 )) 

> min (-H(x, 0) - l E {u £ {x) - ip(x))) > -C. 
xeu 

This shows that ||it e ||ioo is bounded. 

To prove that ||L>ii e ||Loo is bounded independently of e, we first need to prove that H-CH^H^^m 
is bounded by constructing appropriate barriers. 

Let 4> be as in (H2.2). For e small enough, we have that 

$ + H(x, £>$) + 7 £ ($ - ip) < eA$, 

and $ = on dU . Therefore, $ is a sub-solution of (|2.2j) . By the comparison principle, u £ > $ 
in U. 
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Let now d(x) = dist(x, dU). It is well-known that for some 8 > d £ C 2 {Us) and \Dd\ = 1 in 
Ug, where Us :— {x £ 17 : d(ar) < <5}. For /i > large enough, the uniform bound on ||m £ ||£oo 
yields v :— [id > u £ on dUg. Assumption (EGJl) then implies 

v + H(x, Dv) + 7 e (w - ip) - eAv > H(x, /j,Dd) - Cfi > 0, 

for [i is sufficiently large. So the comparison principle gives us that $ < u £ < v in L^. Thus, since 
i/ is the outer unit normal to <9£/, and $ = u £ = w = on dU, we have 

i9w , . du £ . . 9$ , . 

T W < -^-(aO < x , for x e 

w w cw 

Hence, we obtain ||Du e ||£oo( a[/ ) < C. 

\Du £ \ 2 

Next, let us set w £ — — - — . By a direct computation one can see that 

2(1 + (7 £ )')w e + D p H ■ Dw £ + D X H ■ Du £ - (~/ £ )'Du £ ■ Dip = eAw £ - e\D 2 u £ \ 2 . (2.3) 
If ||DM e ||Loo < maxdjDV'llz, 00 , ll-D^IU^sc/)) then we are done. 

Otherwise, max(||£>^||z,Do, ||-Du E ||i,»°(ai7)) < ||£>u £ ||l=. We can choose x 2 £ U such that w £ (x 2 ) = 
max a , g jjw; e (x). Then, using (|2.3[) 

e\D 2 u £ \ 2 {x 2 ) = eAw £ {x 2 ) - 2w £ (x 2 ) - D x H(x 2 , Du £ {x 2 )) ■ Du £ (x 2 ) 

+ ( 7 e )' (Du £ (x 2 ) ■ D^(x 2 ) - \Du £ \ 2 (x 2 )) (2.4) 

< -D x H(x 2 ,Du £ (x 2 )) ■ Du £ (x 2 ). 

Moreover, for e sufficiently small we have 

e\D 2 u £ \ 2 {x 2 ) > e 2 |Aw e (x 2 )| 2 = [u £ {x 2 ) + H{x 2 , Du £ (x 2 )) + j e (u e (x 2 ) - tp(x 2 ))} 2 

(2.5) 

> -\H(x 2 ,Du s (x 2 ))\ 2 -C, 
where we use Lemma l2~2l for the last inequality. Collecting (|2.4p and (|2.5p 

^\H(x 2 , Du £ {x 2 ))\ 2 + D x H(x 2l Du £ (x 2 )) ■ Du £ {x 2 ) < C. 

Recalling hypothesis (EQJl), we must have 

\\Du e \\ L °o = \Du £ (x 2 )\ < C. 

□ 

Thanks to Proposition 12.11 one can show that, up to subsequences, u £ converges uniformly to a 
viscosity solution u of the obstacle problem (|2.ip . 



ADJOINT METHODS FOR OBSTACLE PROBLEMS AND WEAKLY COUPLED SYSTEMS OF PDE 11 



2.2. Proof of Theorem 11.31 We now study the speed of convergence. 
To prove our theorem we need several steps. 

Adjoint method: The formal linearized operator L £ : C 2 (U) — > C(U) corresponding to (|2.2[) 
is given by 

L e z := (1 + (7 £ )> + D p H ■ Dz - eAz. 

We will now introduce the adjoint PDE corresponding to L £ . Let xq g U be fixed. We denote by 
a e the solution of: 



(1 + ( 7 e ) V - div(D p Ha £ ) = sAa e + 6 Xo , 

a e = 0, 



in U, 
on dU, 



(2.6) 



where 8 Xg stands for the Dirac measure concentrated in xq. In order to show existence and 
uniqueness of cr 6 , we have to pass to a further adjoint equation. Let / 6 C'(U) be fixed. Then, we 
denote by v the solution to 



(1 + ( 7 e )> + D p H ■ Dv = eAv + f, 
v = 0, 



in U, 
on dU. 



(2.7) 



When / = 0, by using the Maximum Principle one can show that v = is the unique solution 
to (12. 7p . Thus, by the Fredholm Alternative we infer that (|2.6p admits a unique solution <j e . 
Moreover, one can also prove that a £ € C°°(U \ {^o})- Some additional properties of a e are given 
by the following lemma. 

Lemma 2.3 (Properties of a e ). Let v denote the outer unit normal to dU . Then, 

0o~ £ 

(i) a e > on U . In particular, — — < on dU . 

ov 

(ii) The following equality holds: 



(l + ( 1 £ Y)a £ dx = l + e I -T^dS. 



i)U 



In particular, 







1 

IdU 


dv 



dS < 1. 



Proof. First of all, consider equation (|2.7|) and observe that 

/ > =*> v > 0. 
Indeed, assume / > and let x € U be such that 



(2.8) 



v(x) = minw(x). 
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We can assume that x & U, since otherwise clearly v > 0. Then, for every x € U 

((1 + (Y)'Hx) = eAv(x) + f(x) > 0, 

and ([2T5]) follows, since 1 + (7 s )' > . 

Now, multiply equation f|2 . 6[) by v and integrate by parts, obtaining 



V 



fa £ dx = v(xq). 
Taking into account (|2.8[) . from last relation we infer that 

fa £ dx > for every / > 0, 
and this implies a e > 0. 

To prove (ii), we integrate (|2.6[) over U, to get 

/ (1 + (7 £ )> £ = / div(L> p iJcr e )dx + e / Acr £ dx + 1 

J[7 JU JU 

= [ (D P H ■ v)a E dS + e ( ^dS + l = e [ ^-dS + 1, 
Jdu Jau ov J 9U dv 

where we used the fact that a e = on dU . 

Using the adjoint equation, we have the following new estimate. 
Lemma 2.4. There exists C > 0, independent ofe>0, such that 

1 [ {l + {l e )')\Du £ \ 2 a e dx + e [ |£>V|V dx < C. 

2 Ju Ju 

Proof. Multiplying (|2.3p by a e and integrating by parts, using equation (|2.6[) we get 

1 I (l + h £ Y)\Du £ \ 2 a £ dx + e f \D 2 u e \ 2 a £ dx 

2 Ju Ju 

= -w e (x ) - [ [D X H ■ Du £ - h £ )'DiP ■ Du £ ] a £ dx - e [ w £ ^-dS. 

Ju Jou ov 

Thanks to Lemma |2~51 and Proposition 12. II the conclusion follows. 

Relation (|2.9j) shows that we have a good control of the Hessian D 2 u £ in the support of a 
We finally have the following result, which immediately implies Theorem 11.31 

Lemma 2.5. There exists C > 0, independent of e, such that 
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Proof. Differentiating f|2.2[) w.r.t. e we get 

(2.10) 



(1 + (7 e )'K + d p h ■ Du % +ll = £ A < + A " e , in U. 



In addition, we have u £ (x) — for all x € dU, since u £ (x) = on dU for every e. So, we may 
assume that there exists X2 € U such that |ti|(iE2)| = max x et7 |u|(x)|. 

Consider the adjoint equation (|2.6|) . and choose xq = X2- Multiplying by a £ both sides of (|2.10p 
and integrating by parts, 



Hence, 



By Lemma [2~2l 



itg( x 2) = — / l £ e o- E dx + / Au e a s dx. 
Ju Ju 



K(x 2 )\< ( \j e £ \o- e dx+ [ \Au £ \a £ dx. 
Ju Ju 



u £ — tp , ( u £ — tp 
•^2—7 



(7 £ ) V - V) 



< W 



Hence, thanks to Lemma l2~3l 



while using 



/ |7||cr e dx<C [ (7 £ )V £ da; < C, 

J(7 Ju 



/ f \ 1 / 2 / /" \ V2 

Au £ \a £ dx<(J \Au £ \ 2 a £ dxj I J a £ dx 



C | / |£>V|Vda; 



1/2 



1/2 

cr e da; I < 



C 

rl/2 ' 



Thus, by (l2~TTjl . (I2~13l) and (l2~14)l 



C 



<(a> 2 )| < -772, for s< 1. 



(2.11) 



(2.12) 



(2.13) 



(2.14) 



(2.15) 
□ 



2.3. Dynamic interpretation. We give now a dynamic interpretation of the measure a £ . Thanks 
to the properties given by Lemma [2~31 and arguing as in [EvalOj and |CGT) . we have the following 
theorem. 



Theorem 2.6. There exist 

(i) a measure fi on U x R™, such that 



[ <f)(x,Du £ )cr £ dx^ [_ <f)(x,p)dn, V0G C(UxR n ); 

Ju JuxM n 
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(ii) a measure 71 on U x M. n , and a measure 72 on dU x M. n , such that 



(l + (j E y)cp(x,Du E )(j E dx^ _ 4(x,p)drri, V <p e C(U x 1"), 

U JUxM n 

and 

f r)n£ r 

(x,p)ehf2, y<fre C(dU x R n ), 



-e / ^-cp(x,Du e )dS ^ / 
Vat/ w J dl 



IdUxt 

with the property that 71 (£7 X W n )+-/ 2 (dU xl") = l; 

(iii) a nonnegative definite matrix of measures (m,jk), that we call matrix of dissipation 
measures, such that 

e / 4>{x, Du E )u E x . x .u £ XtXk a E dx -» / cp(x,p) dnijk, V tj> E C(U x E n ), Vj, A; = 1, . . . , n, 

where we used Einstein summation convention; 

(iv) a compact set K C K n such that 

supp^, supp7i, supprTijfe C U x K, supp72 C dU x if. 
Next theorem gives a relation involving the measures fj,, 71, 72, TUkj. 
Theorem 2.7. For any <p £C X (JJ X W 1 ) with <f>(x, •) G C 2 (M") and /or any x eU , we have 
lim 0(x o ,I?ii E (a;o)) = /_ (£> p • (Dip - p) + <p) d 7 i + / cj)d l2 

<^ 0+ JUxR" JdUxWl™ 



Ijk: 



+ /_ ({^' < I > }~ D P ( I ) ' A^>) rf M _ /_ 4>P,Pk dm j 
JUxS." JUxR n 

where the symbol {•, ■} stands for the Poisson bracket, that is 

{F, G} := D P F ■ D X G - D X F ■ D P G, V F,G G C X (U xl"). 

In particular, if <p(x,p) = 4>(x) then 



_ (/)dli+ 4>d-i2+ _ {H, (pjd/j, = (j)(x a ). 

UxR" JdUxR" JUxS." 

Proof. Let us set y> E (x) := (p(x,Du E (x)) for every x G U. Then 

vlt = 4>x t +<>> Pj u € Xi x j > i = l,...,n, 

and 

A(p e = A x (p + 2</> XiPj ul iXj + D p( j> ■ D{Au E ) + (P PjPk u%. x .u%. xlt , 

where we used the notation A x (f> = ^2 i 4>x i x i - Differentiating (|2.2j) w.r.t. x and computing the 
scalar product by D p (p, we get 

(f + (j E y)D p <j> ■ Du E -(Y)'D p <f> ■ Dip + D X H ■ D p <p + H Pj <j> Pi u e XiX . = eD p <P ■ D(Au E ). 
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Thanks to the above calculation on (p x . and Aip £ , 

(1 + ( 7 e )' W ' Du e -(j e )'D p <f> ■ Dip + D X H ■ D p <f> + H Pj - cj> Xj ) 
= e(Atp £ - A x <p - 2<p XiPj u s x . x . - <f> Pm u x . Xj u x . Xk ). 

Hence, adding and subtracting the term D p <p ■ Dip 

(1 + h e Y)D p <P ■ (Du £ -Dip)+D p <p ■ Dip - {H, <P}+D p H ■ Dip 5 
= eAtp £ - £(A x <p + 2(P XiPj ul iXj ) - e4> PjPk u%. x .u%. Xh . 
Multiplying (|2 . 1 6|) by a £ and integrating by parts over U, 



(2.16) 



{D p <p ■ (Du £ -DiP) - tp £ ) (1 + {i £ )')a £ dx+ / (D p( j) ■ Dip - {H, </>}) a £ dx 
u Ju 

- [ tp £ Aw{D p H(j £ )dx~ [ ip £ Aa £ dx+ [ tp £ (l + {~f £ )')<J £ dx (2.17) 
Ju Ju Ju 

= ~ £ f i P E ^~ dS ~ I £ (^x4> + ^4> XlP] u £ x%x P)(J £ dx- e [ <t>p 3Pk u £ x%Xi u £ XiXk (T £ dx. 

JdU ov JU JU 

Recalling equation (|2.6[) and the definition of (p £ , 



(D P 4> ■ {Du £ -DiP) - tp £ ) (1 + (j £ )')a £ dx 
+ [ {D p (p ■ Dip - {H, 4>}) a £ dx + cP(x , Du £ (x )) + £ [ <fi e ^-dS 

Ju JdU OV 

= - / £(A x (p + 2(p XtPj u £ XzX Ja £ dx-£ (p P]Pk u £ XiXi u £ XtXk a £ dx. 
Ju Ju 

Thanks to Holder inequality and using (|2.9p . 

lim / £(A x (P + 2<p x . p .u £ xx .)<7 £ dx <elim / C(l + \D 2 u £ \) o £ dx < lim C-Jl = 0. (2.18) 

s^0+ Ju 1 3 e^O+Jjj e^0+ 

Letting e — > + , using Theorem 12.61 and relation (|2.18p we finally get 



lim (t>(x ,Du s (x ))= / (D P 4> ■ (Dip - p) + <p) d 7l + / 4>d l2 
+ / ({H,(p}-D p <p ■ Dip) dp- / (p PjPk dm jk . 

JuxR n JuxR n 

□ 

2.4. Stochastic process interpretation. In this section we show that the problem could have 
been approached also by using stochastic processes. Let (f2, J 7 , P) be a probability space, and let 
Wt be a n-dimensional Brownian motion on 51. Let e > 0, and let u £ be a solution of (|2.2j) . Let 
T > and consider the solution x e : [0, T] ->■ E" of the SDE 



dx £ = -D p Hhf , Du £ (x e )) dt + V2e dw u 

(2.19) 

x e (0) = x, 
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with x £ U arbitrary. Accordingly, the momentum variable is defined as 

p £ (t) = Du s 

We then define the exit time as 

T§ v := inf{t £ (0, T] : x £ (i) € 017}. 

Let us now recall some basic facts about stochastic calculus. Suppose z : [0,T] 
to the SDE: 



dzi = didt + bij w J t , 



i = 1 



with di and bij bounded and progressively measurable processes. Let ip : R™ x 
function. Then, (p(z,t) satisfies the ltd formula: 



dip = ip Zi dzi + ( (p t + ^bijbjkVziZk ) dt. 



An integrated version of the Ito formula is the Dynkin's formula: 



£[0(z(T))-^(z(O))]=E 
In the present situation, we have 



i n is a solution 



be a smooth 



(2.20) 



OiD«<t>(z(t)) + -bijbjkDi^JMt)) dt 



a t = -D Pt H(x £ ,D x u £ ), b i: j = V2eS i: 



Hence, recalling $TWi and $Z3£§ 

dp £ = D 2 u £ dx £ + eA(Du £ )dt 

= (- D 2 u £ D p H + eA(Du £ ) S j dt + V2eD 2 u £ dw t 

= (d x H + (1 + (7 £ )')p £ - {i £ )'Dip)dt + V2eD 2 u £ dw t , 



where in the last equality we used the identity obtained by differentiating (|2.2p with respect to x. 
Thus, (x e ,p £ ) satisfies the following stochastic version of the Hamiltonian dynamics: 



dx £ = -D p H{x £ ,p £ )dt + V2edw t , 



dp £ = (D x H(x £ ,p £ ) + (1 + (7 £ )')p £ - (Y)'D^ dt + V2eD 2 u £ dw t . 



(2.21) 



We are now going to study the behavior of the solutions u £ of equation (|2.2p along the trajectory 
(x £ (-),p £ (-))- Let cj) £ C^U x R n ). Thanks to the Ito formula and (j2~2T|) . the differential of the 
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function <fi(x. £ (■) , p £ (■)) is given by: 

d<\> = D x 4> • dx. £ + D P 4> ■ dp £ + e (A^ + 2<j> XiPj u e x . x . + 4> PjPk u e XiXj ul %xt ^ dt 
= ({0, H] + (1 + (^)')D p( f> ■ p £ - ( 7 e )'A^ • CV) dt 

+ (sA x </> + 2e(j> xm u e x . x . + e(j> PjPk u%. x .u e x . x ^j dt + V2s (D x <p + D 2 u £ D p (f>) ■ dw t 
= ({0, H} + (1 + h £ )')D P (b ■ (p £ - D^j) + D P 4> ■ D^j dt 

+ (eA x (f> + 2e(f> XiPj u e x . x . + e(j>p jPh ul. x .ul iX A dt + V2e (D x <j> + D 2 u £ D p 4>) ■ dw t 
Thus, by Dynkin's formula we have the following equality 
E[^ £ {T £ du ),p £ {T £ du )) - 4>(x,Du £ {x))] 

" ({0, H} + D p <j) ■ D^j + eA x + 2e4> XiPj u e XiX . + e<t> PjPk u s XiX .u e XiX ^ dt 



E 



E 



1 HI- 



(1 + ( 7 e )' W • (P £ - D^) dt 



which we may also write as 



E 



(D p( j>-(p £ -D4>)-4>)(l + ( 1 £ y)dt 



E 



+ E 



J aU (D p <f>.DTp-{H,ct>} 



dt 



(j)(x ,Du £ (x )) 



0(l + ( 7 £ )')^ 



(2.22) 



E (Tiulp^Tiv))] 



-E 



(eA x 4> + 2e4> XzP3 u £ XiX] + scp p . Ph u e x . x .u e x . Xk ) dt 



Relation (|2.22j) is the analogous of (I2.17[) . and can be as well used (together with suitable estimates) 
to prove Theorem II. 31 

3. Weakly coupled systems of Hamilton-Jacobi equations 

We study now the model of monotone weakly coupled systems of Hamilton-Jacobi equations 
considered by Engler and Lenhart EL91 , and by Ishii and Koike IK91]. For the sake of simplicity, 
we will just focus on the following system of two equations: 

cnMi + C12U2 + Efi(x, Dm) = 

in U, (3.1) 

C21U1 + C22U2 + Ef 2 {x, Du 2 ) = 
with boundary conditions u\ = u 2 = on dU . The general case of more equations or arbitrary 
boundary data can be treated in a similar way. 

We assume that the Hamiltonians Hi,H 2 : U x K™ — > R are smooth satisfying 
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(H3J1) lim I ^\Hj{x,p)\ 2 + D x Hj(x,p) ■ p ) = lim H ^ x ' p ^ = +0O uniformly in x € U, 



|p|-H-oo\2' ) |pK+oo |p| 

for every j = 1,2. Following [EL91] and [IK91] . we suppose further that 
(102) ci2, c 21 <0; 

(E(3J3) there exists a > such that en + Cx2, C21 + C22 > a > 0. 
We observe that, as a consequence, we also have cn,c 22 > 0. Finally, we require that 
(1®4) There exist $1, $ 2 € C 2 (£f) n C X (JJ) with ^ = on dU {j = 1, 2), and such that 

cn$i + ci2$2 + Hi(x, D$i) < in E7, 

c 2 2<i'2+c2i<S>i+ H 2 {x,D<$> 2 ) <0 mU. 
Thanks to these conditions, the Maximum Principle can be applied and existence, comparison and 
uniqueness results hold true, as stated in [EL91] , 

We consider now the following regularized system (here e > 0): 

c\\u\ + C\2U% + H±(x, Du\) = eAu\ 

in U, (3.2) 

c 2 \u\ + C22U2 + H 2 {x, Du e 2 ) = eAu^ 
with boundary conditions u\ = u 2 = on dU . 

Conditions (Ej3]l), (F(3j2) , and (F(3j3) yield the existence and uniqueness of the pair of solutions 

{u\,u%) 



Lemma 3.1. There exists a positive constant C, independent of e, such that 

IKIU-.lluSlk- <c. 

Proof. First of all observe that u\ — u\ — on dU for every e. Thus, it will be sufficient to show 
that u\ and u\ are bounded in the interior of U. 
Let us assume that there exists x € U such that 

max«i(i) = UtCx). 

3=1,2 J 
xeU 



We have 



au\(x) < C\iu\ (x) + ci 2 u e 2 {x) < —Hi(x, 0) < max {—H\{x, 0)) 



where we used (HH3). Analogously, if x € U is such that 

min u e Ax) = uUx), 
3=1,2 J 

then 

u\{x) > — — — u\(x) + — ^ — u e Jx) > - Hl ( x >°> > _ mm(-Hi(x,0)) 

C11+C12 CU+C12 C11+C12 Cn + C12 x&J 
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□ 

Lemma 3.2. There exists a positive constant C , independent of s, such that 

||r>«!l|i»,||i>«|||£»<c. 

Proof. We will argue as in the proof of Proposition l2.il 
Step I: Bound on dU . We shall first show that 

max \Du e Ax)\ < C, 

j=l,2 J 

xedu 

for some constant C independent of e. As it was done in Section^ we are going to construct 
appropriate barriers. For s small enough, assumption (F|2J4) implies that 

cu$i + ci2<I>2 + Hi(x, < eA$i in U, 

c 22 $2 + c 2 i$i + H 2 (x,D<$> 2 ) < eA$ 2 in U, 

and $1 = $2 = on dU. Therefore, (i>i,$ 2 ) is a sub-solution of (|3.2j) . By the comparison 
principle, u £ > <£>., in U , j — 1, 2. 

Let 5, and Us be as in the proof of Proposition 12. II For // > large enough, the uniform 

bounds on ||uf ||l°o and HwllU 00 yield v :~ /id> u £ on dUg, j = 1,2, so that 

(en + Ci2)v + H\(x, Dv) — sAv > Hi(x, [iDd) — [iC in U, 
(c 2 i + Ci2)v + H 2 (x, Dv) — sAv > H 2 (x, fxDd) — fiC in U. 
Now, we have <F,- = u £ = v = on <9[7. Also, thanks to assumption (E(311), for n > large enough 

(en +ci 2 )« + fl'i(x, Dd) -eAw > mU, 

(c 21 +c 22 )v + H 2 (x,Dv)~eAv>0 in U, 

that is, the pair (v,v) is a super-solution for the system (|3.2p . Thus, the comparison principle 
gives us that <&j < u| < « in Us, j = 1, 2. Then, from the fact that $ 3 = = u = on <9[/ we get 

^W^W^W, for xe dU, j = 1,2. 
cw w w 

Hence, we obtain ||-Dw|||i,°o(3[/) < Cj J ' = 1)2- 

Step II: Bound on {/. 

\Du £ \ 2 

Setting us? = — ^ — , j = 1, 2, by a direct computation we have that 

2cnwf + DpHt ■ Dwl + c 12 Du\ ■ Du e 2 + D X H X ■ Du\ = eAw\ - e\D 2 u\ | 2 , 

(3.3) 

2c 22 W2 + D P H 2 ■ Dw\ + c 2l Du\ ■ Du% + D X H 2 ■ Du e 2 = eAw £ 2 - e\D 2 u £ 2 \ 2 . 
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Assume now that there exists x £ U such that 

maxwf(i) = wUx). 

3 = 1,2 J 

xeu 

Then, we have 

e\D 2 u\\ 2 {x) = eAwf(x) ~ 2c u «jf (x) - c 12 Du\{x) ■ Du £ 2 (x) - D X H X ■ Du\{x) 
< -2(c u + c 12 )w\{x) - D x Hi ■ Du\{x) < -D x Hi ■ Du\{x). 
Now, for e sufficiently small 

e\D 2 u\{x)\ 2 >e 2 \Au\{x)\ 2 = [c u u\ {x) + c 12 u%{x) + Hx{x,Du\ (x))] 2 
^^H^DuKx^-C. 
Collecting the last two relations we have 

i|Hi(2,£>uf (x))| 2 + D x H 1 (x,Du{(x)) ■ Du\{x) < C. 

Recalling condition (Ej3]l) the conclusion follows. □ 



Adjoint method. At this point, we introduce the adjoint of the linearization of system (|3.2 
The linearized operator corresponding to (13.21) is 



L 6 (z 1 ,z 2 ) := 



D p H 1 (x, Du\) ■ Dzi + cnzi + c 12 z 2 — eAzi, 
D p H 2 (x,Du 2 ) ■ Dz 2 + c 22 z 2 + c 2 izi - eAz 2 . 
Let us now identify the adjoint operator (L e )* . For every u 1 , v 2 £ C^°(U) we have 

((L E )>V 2 ), (zu**)) ■■= ((" V), L% Zl ,z 2 )) 
= {v\[L s {z u z 2 )} 1 ) + {v 2 ,[L^z 1 ,z 2 )] 2 ) 

= / [D p Hi(x, Dul) ' Dz\ + c\\Z\ + c\ 2 z 2 — sAz\\ v 1 dx 
Ju 

+ / [D p H 2 (x, Dul) ' D z z + c 22^2 + c 2 izi - sAz 2 ] v 2 dx 
Ju 

= / [-div(D p i?i^ 1 ) + cuv 1 + c 21 v 2 - eAv 1 ] z x dx 
Ju 

+ / [-~&iv{D p H 2 v 2 ) + c, 22 v 2 + c^v 1 - eAv 2 ] z 2 dx. 
Ju 

Then, the adjoint equations are: 

-div(L>piJicj 1 ' e ) + oner 1 ' 6 + c 21 a 2 ' £ = eAa 1 ' 6 + (2 - i)8 Xo in U, 
~div{D p H 2 a 2 ' £ ) + c 22 a 2 ' £ + c 12 a he = eAa 2 ' e + (i - l)S Xo in U, 



(3.4) 
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with boundary conditions 

on dU, 
on dU, 

where i £ {1, 2} and xo € U will be chosen later. Existence and uniqueness of a 1 ' 6 and <7 2 ' e follow 
by Fredholm alternative, by arguing as in Section [2j and we have a 1,£ ,a 2,£ € C°°(U \ {xq}). We 
study now further properties of cr l£ and a 2 ' £ . 




Lemma 3.3 (Properties of a 1,£ , a 2,£ ). Let v be the outer unit normal to dU. Th 

— da j ' e 

(i) a 3,6 > on U. In particular, — — < on dU (j — 1.2). 

ov 

(ii) The following equality holds: 
2 ,. ,. Qgj 



j ■ 

In particular, 



E ( + c *y ,£ dx - e L ~&r dS ) = 1 - 



U JdU 



y~] / {cji + c j2 )o- j ' E dx < 1. 



(3.5) 



Proof. First of all, we consider the adjoint of equation (|3.4[) : 

D p Hi(x,Du € 1 ) ■ Dzi + cnzi + C12Z2 - eAzi = / 1; 
D p H 2 {x, Du e 2 ) ■ Dz 2 + C22Z2 + C21Z1 - eAz 2 = /2, 
where /i,/a € C(U), with boundary conditions zi = z 2 = on 9f7. Note that 

/i,/ 2 >0=» rnin ^(x) > 0. (3.6) 

Indeed, if the minimum is achieved for some x € dU, then clearly Zi, £2 > 0. Otherwise, assume 

min z-j(x) = Zi(x), 
3=1,2 J 

for some x € U. Using condition (Ej3j2) 

(en + cu)zx(x) > cnzi(x) + c 12 z 2 {x) = eAzi(x) + fi{x) > 0. 

Thanks to (IEJ3), (EH) follows. 

Let us now multiply (|3.4[) i and (|3.4I) 9 by the solutions z\ and z 2 of (|3.5p . Adding up the 
relations obtained we have 

f 1 a 1 ' £ dx + / f 2 a 2 ' £ dx = (2 - i)z 1 (x ) + (1 - i)z 2 (x ). 
u Ju 

Thanks to (|3.6[) . from last relation we conclude that 

f x o x ' £ dx+ f f 2 a 2£ dx > 0, for every f u f 2 > 0, 
u Ju 
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and this implies that a 1 ' 6 \a 2 - e > 0. To prove (ii), it is sufficient to integrate equations p.4ft i and 
(I3.4[) 9 over U, and to add up the two relations obtained. □ 

Lemma 3.4. There exists a constant C > 0, independent of e, such that 

el |D a vf|V' B dx + e [ \D 2 u £ 2 \ 2 cr 2 > £ dx < C. 
Ju Ju 

Proof. To show the lemma, one has first to multiply equations (|3.3[) i and (|3.3|) q by a 1 ' 6 and a 2 ' e 
respectively. Then, adding up the relations obtained and using (|3.4[) , thanks to Lemma 13.11 and 
Lemma l3~2l the conclusion follows. □ 

We now give the last lemma needed to estimate the speed of convergence. Here we use the 
notation u e £ (x) :— duj(x)/de, j = 1,2. 

Lemma 3.5. There exists a constant C > 0, independent of e, such that 

max|4 £ (x)| < 

xeu 

Proof. Differentiating (|3.2I) w.r.t e we obtain the system 

ci\u\ e + ci 2 u 6 2 e + D p Hi ■ Du\ e = eAu\ e + Au{ , 

(3.7) 

C2\u\ e + c 2 2"2,£ + D p H 2 ■ Du e 2 e = sAu E 2 e + Au 2 . 
Since u\ £ — u 2 s — on dU, we have 

max ii, Jx) — max u% Ax) = 0. 

xedu 6 xedu ' e 

Assume now that there exists x £ U such that 

m ax|u| (x)\ = \uf (x)\, 
xeu 

and let a 1 ' 6 , a 2,E be the solutions of system Q3.4p with i = 1 and xq = x. 

Multiplying equations f|3.T[) i and (|3. 7^ 9 by cr 1:£ and a 2 ' e respectively and adding up, thanks to 
(13.4)) we obtain 

u i,e(%) = / Au\a 1 ' £ dx+ I Au E 2 a 2 ' £ dx. 
Ju Ju 

Thanks to Lemma 13.41 and repeating the chain of inequalities in (|2.14l) one can show that 

C 



Auj a°' e dx 
u 



and from this the conclusion follows. □ 
We can now prove the following result on the speed of convergence. 
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Theorem 3.6. There exists C > 0, independent of e, such that 

IK - «i|U,°°, \\u e 2 - u 2 \\ L ™ < Ce 1/2 . 
Proof. The theorem is a direct consequence of Lemma 13.51 □ 

4. Cell problem for Weakly coupled system of Hamilton-Jacobi equations 

In this section, we study the following weakly coupled systems of Hamilton-Jacobi equations: 
c\U\ — c\u 2 + Hi(x, Dui) = Hi 

in IT, Hi,H 2 eR, (4.1) 

-c 2 ui + c 2 u 2 + H 2 (x, Du 2 ) = H 2 
which is the analog of the cell problem for single equation introduced by Lions, Papanicolaou, and 
Varadhan [LPV88] . We will assume that H U H 2 G C°°(T" x E n ), and: 

lim [ ~\Hj(x,p)\ 2 + D x Hj(x,p) -p- 8nc 2 \p\ 2 ) = +oo uniformly in x G T™, j = 1,2; 
(H32) ci.ca >0. 

It is easy to see that the coefficients of u\ , u 2 in this system do not satisfy the coupling assumptions 
of the previous section. Indeed, as it happens for the cell problem in the context of weak KAM 
theory, there is no hope of a uniqueness result for (|4.1[) . 

To find the effective Hamiltonians Hi,H 2 we use the same arguments as in |Traj . First, for 
every s > 0, let us consider the following regularized system: 

(ci +e)u\ ~ c x u\ + Hi(x,Dul) = e 2 AuJ 

in TP. (4.2) 

(c 2 + e)u 6 2 - c 2 u\ + H 2 (x, Du 2 ) = e 2 Au 2 : 

For every e > fixed, the coefficients of this new system satisfy the coupling assumptions (102) 
and (Edl3) of the previous section. Thus, (|4.2[) admits a unique pair of smooth solutions u\,u\. 
In particular, this implies that u\ and u 2 are T™-periodic. 
The following result gives some a priori estimates. 

Theorem 4.1. There exists C > 0, independent of e, such that 

||eu!||L-,||eu|||i-,||I>«!l|£-.,||I>t4||i- < C. 

Proof. Our proof is based on the Maximum Principle. Without loss of generality, we may assume 
that 

max{ew<;(ir)} = eu\ (x E ). 
for some Xq G T". Applying the Maximum Principle to the first equation of (14.21) . 



eu\(xl) < (ci + e)uf(x £ ) ~ d^(xg) < ~H'(x s Q ,0) < C, (4.3) 
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and this shows the existence of a bound from above for eu\ and eu 2 . Using a similar argument 
one can show that there is also a bound from below, so that 

||eu;||£~,||eu!|| £ « <C. (4.4) 

We observe that the previous inequality doesn't provide any bound for the difference u\ (x) — u 2 (x) 
in a generic point x G T". Nevertheless, thanks to (|4.4[) we have 

u\{xl)-ul{x%)<--H 1 {x%,{))--ul{xl)<C. (4.5) 

Ci Ci 

Then, (g3]) and (g3]) imply that 

Cl \ul{xl)-u%[xl)\<C. (4.6) 
\Du e ,\ 2 

In order to find a bound for the gradients, let us set w 6 j — — ^ — , j — 1,2. Then, by a direct 
computation one can see that 



2(ci + e)w{ + D P H 1 ■ Dw\ - c x Du\ ■ Du\ + D X H X ■ Du\ = e 2 Aw{ - e 2 \D 2 u\\ 2 
2(c 2 + e)w e 2 + D p H 2 ■ Dw e 2 - c 2 Du\ ■ Du e 2 + D X H 2 ■ Du e 2 = e 2 Aw e 2 - e 2 \D 2 u e 2 \ 2 
Without loss of generality, we may assume that there exists x\ G T™ such that 

max \w e Ax)\ — w s (x e ). 

Then, by the Maximum Principle 

e 2 \D 2 u\{x\)\ 2 < -2(ci + e)w\{xl ) + Cl Du\{x\) ■ Du e 2 {x\) - D X H X ■ Du\(x\) 
< -DxH^DuKxl). 
Moreover, for e sufficiently small 



m 



T" 



(4.7) 



e 2 \D 2 u\{x\)\ 2 > e\Au\{x\)) 2 = [H^xl , Du\ (x\ )) + ( Cl + e)u\(x\) - Cl u%{x\)f . (4.8) 
Also, thanks to (O)) and (HI?)) 

\{ci + e)u\{x\) - ciu e 2 (xl)\ 

< e\ul{x\)\ + Cl \u\{x\) - u\{xl)\ + Cl \ul{x\) - u e 2 {x e )\ + ci|uf (ig) - u E 2 {x s )\ 
<C + c x \u\ (xf ) - u\{xl)\ + ci|u|(xf) - u|(as§)| 

< C + 2ci\Dul(xl)\\xl - x E \ < C + 2ciV^\Dul(xl)\, 
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where we used the fact that the diameter of T" is y/n,. Last relation, together with (14.81) . gives 
that 

e 2 \D 2 u\{x\)\ 2 > [Hxixl, Du\{xD) + (ci + e)u\{x\) - c^xl)} 2 
> ^\Hi{x\ ,Dul(xf))\ 2 -C- 8nc 2 |i^(^)| 2 . 

Using last inequality and (|4.7[) we have 

^H^xuDuKxl))^ + D X H X ■ Bu\(x\) - Snc 2 \Du\{x\)\ 2 < C. 

This, thanks to condition (bEll): gives the conclusion. □ 

Thanks to Theorem 14. 11 up to subsequences, 

eu\ — > Hi, eu 2 — > H 2 , uniformly in T n , 

for some constants if 1, if 2 S K- Furthermore, still up to subsequences, 

uf — mintt? — > Ui, u% — mintto — > u 2 , uniformly in T™, 

where U\ and u 2 are viscosity solutions of (|4.ip . 

In general, Hi and H 2 are not unique. Indeed, let ui,u 2 be viscosity solutions of (I4.1[) . Then, 
for every pair of constants Ci, C 2 , the functions ui := ui + Ci and u 2 := u 2 + C2 are still viscosity 
solutions of (|4.ip . with new effective Hamiltonians 

Hi=Hi+ Cl (Ci - C 2 ), H 2 =H 2 + c 2 {C 2 - Ci). 

Anyway, we have C2H1 + ciH 2 = C2H1 + ciH 2 . This suggests that, although Hi and H 2 may 
vary, the expression c 2 Hi + ciH 2 is unique. Next theorem shows that this is the case. 

Theorem 4.2. There exists a constant (jgl such that 

c 2 Hi + ciH 2 = /1, 

for every pair (Hi,H 2 ) G R 2 such that the system (|4.1[) admits viscosity solutions ui,U2- 

Proof. Without loss of generality, we may assume ci = c 2 = 1. Suppose, by contradiction, that 
there exist two pairs (Ai,A2) S R 2 and (ni,^) € K 2 , and four functions ui,U2,ui,U2 € C(T") 
such that Ai + A2 < /ii + /i2 and 

— u 2 + Hi (x, Dui) = Ai 

inT'\ 

—til + u 2 + H 2 (x, Du 2 ) = X 2 
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and 

ui - u 2 + Hi(x,Dui) = fii 

in T™. 

-Ml +U2 + H 2 (x,Du 2 ) = 

By possibly substituting u\ and M2 with functions u\ := Ux + C\ and U2 := 1*2 + C2, for suitable 
constants C\ and C2, we may always assume that Ai < fj,i, A2 < \i 2 . 

In the same way, by a further substitution u\ := ui + C3, U2 := «2 + C3, with C3 > large 
enough, we may assume that u\ > ui, U2 > u 2 . Then, there exists e > small enough such that 

(e + \)u\ - u 2 + Hi(x, Dui) < (e + l)ui ~u 2 + H\(x,Du\) 

in T". 

(e + l)u 2 - ui + H 2 (x,Du 2 ) < (e + l)u 2 — ui + H 2 (x,Du 2 ) 

Observe that the coefficients of the last system satisfy the coupling assumptions (E02) and (tJ3]3). 
Hence, applying the comparison theorem in [EL91] and |IK91j . we conclude that u\ < u\ and 
u 2 <u 2 , which gives a contradiction. □ 

In the sequel, all the functions will be regarded as functions defined in the whole R™ and 
Z™-periodic. Next lemma provides some a priori bounds on u\ and til- 
Lemma 4.3. There exists a constant C > 0, independent of s, such that 

K(x) - u\{y)\, \ul{x) - ul{y)\, \u\{x) - u £ 2 (y)\ <C, x, y e M". 

Proof. The first two inequalities follow from the periodicity of u\ and u 2 , and from the fact that 
Du\ and Du 2 are bounded. 

Let us now show the last inequality. Without loss of generality, we may assume that there 
exists xo € T n such that 

max {ufa)} = u\{xl). 

x£T" 

Combining the second inequality of the lemma with (14.51) . 

ul{x) - u s 2 (y) < ul(x e ) - u e 2 (x e ) + u e 2 (x e ) - u s 2 (y) < C, x,yeM. n . 
The proof can be concluded by repeating the same argument for min^i^ {u s Ax)} . □ 
The following is the main theorem of the section. See also [Tra for similar results. 

Theorem 4.4. There exists a constant C > 0, independent of e, such that 

Weul-H^Loo^eul-H^ <Ce. 
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Adjoint method: Also in this case, we introduce the adjoint equations associated to the 
linearization of the original problem. We look for er 1,£ , <7 2,e which are T ra -periodic and such that 

-dxv(D p H 1 tr 1 ' e ) + (ci + e)<r l ' e - c 2 a 2 ' £ = e 2 Ao- 1 ' e + s(2 - i)5 Xo in TP, 



(4.9) 

-div(D p H 2 a 2 ' s ) + (c 2 + e)a 2 ' £ - c x a x < e = e 2 Aa 2 - £ + s(i - l)5 Xo in T" 



where i € {1,2} and xo € T n will be chosen later. The argument used in Section [5] gives also in this 
case existence and uniqueness for a 1,£ and a 2,£ . As before, we also have a 1 ' 5 , u 2 - £ £ C°°(T n \{xo}). 
The next two lemmas can be proven by using the same ways as in Lemma 13.31 and Lemma [ 



Lemma 4.5 (Properties of a 1,£ ,a 2,£ ). The functions a ' e ,a ' s satisfy the following: 

(i) a^ £ > on T" (j = 1,2); 

(ii) Moreover, the following equality holds: 



Lemma 4.6. There exists a constant C > 0, independent of s, such that 



e 2 



] f |U 2 «;|V' e dx < C, 

e 2 J \D 2 u £ \ 2 a 2 - £ dx < C. 
Finally, next lemma allows us to prove Theorem 14.41 
Lemma 4.7. There exists a constant C > 0, independent of e, such that 

max|(euf) s |, rnax|(eu|) E | < C. 

Proof. Differentiating (|4.2I) w.r.t. e, 

D p Hi ■ Du £ l e + (ci + e)u\ E + u\- ciu £ 2 e = e 2 Au £ l £ + 2eAu\, 

D P H 2 ■ Du £ 2 e + (c 2 + e)u £ 2 e + u £ 2 - c 2 u\ E = e 2 Au £ 2 e + 2eAu £ 2 , 
where we set e := du^/de, j = 1, 2. Without loss of generality, we may assume that there exists 
x 2 £ T™ such that 

max(eu £ (x)) £ — max {su £ £ (x) + u £ (x)} = eu\ e (x 2 ) + u\{x 2 ). 

Choosing xq — x 2 in the adjoint equation (|4.9p . and repeating the same steps as in Theorem 12.51 
we get 



eu £ l £ (x2) + j u\<j 1,e dx+ I u 2 o- 2 ' £ dx 
<2e j lAuWo- 1 ^ dx + 2e [ \Au £ 2 \cr 2 ' £ dx < C, 



(4.10) 
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where the latter inequality follows by repeating the chain of inequalities in (|2.14j) and thanks to 
Lemma 14.61 Using Lemma 14.31 and property (ii) of Lemma 14.51 we have 



u\(x)<j ' 6 dx + / u £ 2 {x)<t ,e dx — u\{x 2 ) 



{u\(x) — (x 2 ))cr 1 ' e dx + / (u e 2 (x) — u\{x 2 ))a 2 ' £ dx 



< C. 



In view of the previous inequality, (14. 10)) becomes 



eu\ E {x 2 ) +u{(x 2 ) < C, 



thus giving the bound from above. The same argument, applied to min,, = i 2 {£u £ {x)) e , allows to 
prove the bound from below. □ 



Proof of Theorem \4-4\ The theorem immediately follows by using Lemma 14.71 



□ 



5. WEAKLY COUPLED SYSTEMS OF OBSTACLE TYPE 

In this last section we apply the Adjoint Method to weakly coupled systems of obstacle type. 
Let Hi,H 2 : U x R™ — s- R be smooth Hamiltonians, and let ipi,^ '■ U —> R be smooth functions 
describing the obstacles. We assume that there exists a > such that 



tpi, ip2>a in U, 



(5.1) 



(5.2) 



and consider the system 

max{ui — u 2 — ipi,Ui + Hi(x,Dui)} = in U, 

max{u 2 — u\ — ip2, u 2 + H 2 (x, Du 2 )} = in U, 

with boundary conditions u\ \du= u 2 \du= 0. We observe that (|5.1|) guarantees the compatibility 
of the boundary conditions, since ^p\,ip 2 > on dU . 

Although the two equations in (|5 .3[) are coupled just through the difference u\ — u 2 , this problem 
turns out to be more difficult that the correspondent scalar equation (|2.ip studied in Section [5J 
For this reason, the hypotheses we require now are stronger. We assume that 

(£01) Hj(x, ■) is convex for every x € U, j = 1, 2. 
(L^5j2) Superlinearity in p: 

Hj{x,p) 



lim 



oo uniformly in x, j = 1, 2. 



M-s-oo \p\ 

(H33) \D x H 3 {x,p)\ < C for each (x,p) eUx R™, j = 1,2 
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(H34) There exist $ x , $ 2 € C 2 (C/) n C^U) with = on dU (j = 1, 2), -r/> 2 < $i - $2 < V'l, 
and such that 

$j + £>$,•) < in U (j = 1, 2). 

Let e > and let 7 e : K — >• [0, +oo) be the function defined by (|1.7[) . We make in this section the 
additional assumption that 7 is convex. We approximate (|5.3[) by the following system 

u\ + H 1 (x 7 Du £ 1 ) + j £ (u\ - 1*2 - «/>i) = eAu^ in J7, 

+ H 2 (x,Du £ ,) + 7 e (it| - uf - fa) = eAuf in U. 

We are now ready to state the main result of the section. 



(5.3) 



L £ (z u z 2 ) 



Theorem 5.1. There exists a positive constant C, independent of e, such that 

max{||uf - UiHio, ||w| ~ u 2 \\l™} < Ce 1/2 . 
In order to prove the theorem we need several lemmas. In the sequel, we shall use the notation 
6\ := u\ — 1*2 — ipi , 6*2 := u 2 — u\ — tp 2 - 
The linearized operator corresponding to (|5.3p is 

zi + D v H 1 (x,Du\) ■ Dzi + (7 £ )' |o; {zi - z 2 ) - eAzi, 
z 2 + D p H 2 (x,Du e 2 ) ■ Dz 2 + ( 7 £ )' | . (z 2 - zi) ~ eAz 2 . 
Then, the adjoint equations are: 

(1 + ( 7 £ )' \ei)cr l ' e - dWiDpH.a 1 ^) - ( 7 £ )' \e S <r 2 > £ = eAa 1 ^ + (2 - i)5 Xo in U, 

(5.4) 

(1 + ( 7 £ )' \eg)o- 2 > e - dW(D p H 2 a 2 ' £ ) - ( 7 £ )' | 0f a 1 ^ = £ Aa 2 < £ + (i - l)6 Xo in (7, 
with boundary conditions 

on dU, 
on 9/7, 

where i € {1,2} and xq E U will be chosen later. By repeating what was done in Section [51 we 
get the existence and uniqueness of a x,e and a 2 ' e by Fredholm alternative. Furthermore, a 1,£ and 
a 2,£ are well defined and a 1,£ ,a 2,£ € C°°(U \ {^o})- In order to derive further properties of a 1 ' 5 
and cr 2 ' e , we need the following useful formulas. 




Lemma 5.2. For every tp\,ip 2 € C 2 (U) we have 

- f r (l £ )'k Via 2 - e dx 

(5.5) 



(2 - i)^ (x ) = -e [ ^— dS- [ ( 7 £ )' | 9 . ^ a 2 ' 6 
Jot/ w Jrj 



+ / [(1 + (7 £ )' Vi + DpHi ■ D Vl - eA Vx ] a 1 ' 6 dx, 
Ju 
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(i - l)ip 2 (x ) = -e 



da 



2.s 



av 



du 



<p 2 dS- ( 7 e )' |o; V2<?' e dx 



+ I [(1 + (7 £ )' \o{)<P2 + D P H 2 ■ Dcp 2 - eA V2 ] a 2 ' 6 dx, 
Ju 



where v is the outer unit normal to dU . 



Proof. The conclusion follows by simply multiplying by tpj (j — 1,2) the two equations 
and integrating by parts. 



We can now prove the analogous of Lemma 12.3 



Lemma 5.3 (Properties of a > e ,a > e ). Let v be the outer unit normal to dU. Then 

— da je 
(i) G 2 ' e > on U . In particular, — — < on dU (j = 1,2). 



(ii) The following equality holds: 



E 

3=1 



In particular, 



a J ' dx — e I — — 

u Jdu dv 



2 r 

J ^' e dx < 1. 



dS 



Proof. First of all, we consider the adjoint of equation (J57 

2i + D p Hi{x, Du\) ■ Dzi + ( 7 e )' | e e (zi - z 2 ) - eA Zl = fx, 
z 2 + D p H 2 {x, Du e 2 ) ■ Dz 2 + ( 7 e )' \ge (z 2 - z x ) - eAz 2 = f 2 , 
where fx, f 2 € C(U), with boundary conditions z\ = z 2 = on dU . Note that 



/i, h > => min zj(x) > 0. 



Indeed, assume 



min z,-(x) = zi(x), 

3 = 1,2 J 
x£U 

for some x € U. Then, 

Zl {x) > zx(x) + ( 7 e )' 1 i (zi(x) - z 2 {x)) = eA Zl (x) + fx(x) > 0. 
Adding up relations (|5.5p and (|5.6[) with y?i = zx and = z 2 we get 



f 1( j L - e dx+ / f 2 a z - £ dx = (2-i)zx{x ) + (1 -i)z 2 (x ). 

(7 J(7 



Thanks to (15.81). from last relation we conclude that 



/i,/ 2 >0 



/io- 1 ' 5 + / f 2 a^ e dx > 0, 
u Ju 
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and this implies that c 1 ^, a 2 - e > 0. 

The proof of property (ii) follows by choosing tpx = (f2 = 1 in (|5.5[) and (|5.6p . and summing 
the relations obtained. 

We are now able to prove a uniform bound on u\ and tt|. 
Lemma 5.4. There exists a positive constant C , independent of e, such that 



KlU-JuIIU" <c. 



Proof. If the maximum xsxaxj—i l 2{Uj(x)} is attained on the boundary, then 



xeu 



maxlu e Ax)\ = maxlv5(x)} = 0. 

j=l,2 J 3=1,2 ■> 

xeu xedu 



Otherwise, assume that there exists x € U be such that 



max! u e Jx 
xeu 



Then, 



u\(x) = eAu\{x) - Ht{x, 0) - i £ {u\(x) - u%{x) - ipi{x)) 



< -H x {x,Q) < sup (-ifi(ar,0)) < C. 



xeu 



Analogously, to prove the lower bound suppose that 



min \u e Ax 

xeu 




for some x € U. Then, u\{x) — u E 2 {x) — ipi{x) < and so 



ul(x) = eAul(x) - Hi(x,0) - i £ {u\{x) - u e 2 (x) - ^i(x)) 



eAu\(x) - #i(x,0) > -Fi(x,0) > inf (-i?i(x,0)) > C*. 



Next lemma will be used to give a uniform bound for Du\ and Dm|. 



Lemma 5.5. We have 




where C is a positive constant independent of e. 
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Proof. First of all, observe that condition (L(5jl) implies that 

H j (x,p)-D p H j (x,p)-p<H j (x,0), for every (x,p) <E~UxW\ j = 1,2. 
In the same way, the convexity of 7 implies 

7 £ (*) " [(7 £ )'(*)] s = 7 (J) ~ [(7) (£)] ~ < 7(0) = 0. 

Equation (|5.3p i gives 

= uf + ifi(a;,Dui) + 7 e | e i -eAuf 
= uf + D p H 1 (x, Du\) ■ Du\ - eAu\ + H x (x, Du\) - D p Hx{x, Du\) ■ Du\ 

+ 7 £ k "(7 £ )' le| ^ + (7 e )' \ei K - u\) - ( 7 £ )' \ el V'l- 
Multiplying last relation by cr 1 ' 6 , integrating and using (|5.9I) and (|5.10|) 

/ ( 7 e )' | e i Vicr 1,£ d2;= / [H\{x, Du\) — D p H 1 (x, Du\) ■ Du\] a 1 ' 6 dx 
Ju Ju 

+ [ [7 £ k-(7 £ )' Wjl] ^ £ dx 
Ju 

+ I [(1 + (7 £ )' leO u\ + D p H x {x, Du\) ■ Du\ - eAu\ - ( 7 e )' \ g i u§] a l e dx 
Ju 

< [ H 1 (x,0)a 1 ' £ dx 
Ju 

+ /[(! + (7 £ )' lei) u\ + DpHtix, Du\) ■ Du\ - eAu\ - ( 7 e )' U u £ 2 ] a l e dx. 
Ju 

Analogously, 

f ( 7 £ )' \ 9 2^ 2 <T 2 ^dx < [ H 2 (x,0)a 2 ' e dx 
Ju " Ju 

+ I [(1 + (7 £ )' |e|) u\ + D p H 2 ( Xl Dul) ■ Du\ - sAu% - ( 7 £ )' | ; u|] a 2 > £ dx 
Ju 

Summing up the last two relations and using (|5.5[) and (|5.6[) 

( 7 e )' | fl i ^ dx + [ ( 7 e )' \ 8 2 i> 2 a 2 - e dx<{2- i)ul(xo) + (i - l)«I(a:o) 

+ ||ffi (-,0)11^00 f o-^ e dx+\\H 2 {;0)\\ L ^ / <r 2 ' e cfe. 

7 [7 Ju 



Thus, 



( 7 e )' lei (T^darH- J ( 7 e )' 1 0i CT 2 ^dx< ^— -u\{x )+ l —^ul{xo) 



a ju 

where we used (|5.ip . Lemma [5751 and Lemma [57^1 

We can finally show the existence of a uniform bound for the gradients of u\ and u 2 . 
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Lemma 5.6. There exists a positive constant C, independent of e, such that 

puflU-, \\Dul\\ L ~<C. 

Proof. Step I: Bound on dU. 

As it was done in Section[3J we are going to construct appropriate barriers. For e small enough, 
assumption (EJ3]4) implies that 

$1 + Hi^-D^i) +7 £ ($i - $2 - ipi) < eA$i in U, 
$ 2 + H 2 (x, L>$ 2 ) + 7 e ($2 - $i - V>2) < eA$ 2 in C/, 

and $i = $2 = on dU. Therefore, ($i,$2) is a sub-solution of (|5.3[) . By the comparison 
principle, Uj > $j in [/, j = 1, 2. 

Let 8, and C/5 be as in the proof of Proposition 12. II For /1 > large enough, the uniform 

bounds of ||uf||it» and HulIU 00 yield v := fid > u £ j on dUs, j = 1,2, so that 

v + H 1 (x 7 Dv) +7 £ (v -v-ijjx)- eAv = v + Hi{x,Dv) - eAv > H 1 (x 7 ^Dd) - \iC in U, 
v + H 2 {x 1 Dv) +7 £ (w -v-ip 2 )- sAv = v + H 2 {x,Dv) - eAv > H 2 {x 1 ^Dd) - \iC in U. 

Now, we have <£>j = u £ = v = on dU. Also, thanks to assumption (HSJ2), for \i > large enough 

{v + Hi(x,Dv) +7 e (-u -v-ipi)- eAv > in U, 
v + H 2 (x,Dv) +7 e (w -V-1P2)- eAv > in U, 

that is, the pair (v,v) is a super-solution for the system (|5.3I) . Thus, the comparison principle 
gives us that $j < < i>j in Us- Then, from the fact that &j = = v = on dU we get 

^(x) < -^x < ^(x), for z€ 3*7. 

CW CW 

Hence, we obtain ||-DUj||L°°(3[/) <C, j = 1, 2. 
Step II: Bound on U. 

Assume now that there exists x <E U such that 

max w j (x) — wl(x), where Wj(x) := ^\Duj\ 2 , j — 1,2. 

By a direct computation one can see that 

2(1 + (7 £ )' IfljM + A>#i • f t«f + D X H X ■ Du\ - ( 7 £ )' \ 9 . Du\ ■ (£% + Du £ 2 ) = eAw\ - e\D 2 u\ | 2 . 
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Multiplying last relation by a 1 ' 6 and integrating over U 

2 / w\a 1 > e dx + / D V H X ■ Dwlcr 1 ' 6 dx - e / Awf a 1 ' 5 dx + / e\D 2 u\ |V< £ dx 
Ju Ju Ju Ju 

+ / D X H X ■ Du\o x * dx + \ f ( 7 e )'k [\Du\\ 2 + \Du\- Du' 2 \ 2 -\Du e 2 \ 2 } a 1 ' 6 dx 
Ju * Ju 

- I il e )' k Du\-Di> 1 (r l ' e dx = Q. 
Ju 

Then, using equation (I5.5[) i with i — 1 and xq = x 

[ w £ 1 a 1 - £ dx+ [ e\D 2 u\\ 2 a^ e dx - [ ( 7 e )' | . Du\ ■ D^a 1 " dx 
Ju Ju Ju 

+ { D X H X - Du\a l > e dx + \ [ ( 7 £ )' \ 8{ [\Du\ - Du e 2 \ 2 - \Du £ 2 \ 2 } a 1 ' 5 dx 
Ju - J i 



(5.11) 



(5.12) 



w\{x) + e / -^w\ dS + / ( 7 e )' | . w\ a 2 ' £ dx = 0, 
Jdu v v Ju 



which implies 

wf(x)- [ ( 7 £ )' \e t w e 2 a 1 ' s dx+ f ( 7 £ )' | 9 « wf a 2 ^ dx 
Ju Ju 

< [ (7 £ )' k Du\ ■ Difjia 1 ' 5 dx - [ B X H X ■ Du\a l ' e dx - e [ ^—wfdS. 

Ju Ju JdU VV 

Let now rj > be a constant to be chosen later. Using Step I and Lemmas 15.31 and 15.51 thanks to 
Young's inequality 

/ (7 £ )' k <a 1 ' £ dx+ /( 7 £ )' | fl . w\a 2 ^dx 
Ju Ju 



< / (7 E )' k 
Ju 



\\DH\h 



2if 



a 1 ' 6 dx 



\D x Hi\\ L o 



2 if 



<tf(C + l)wl(x) + C[ 1 + ^ 



(5.13) 



a^dx + C 



In the same way, considering the analogous of equation (|5.1ip for the function w 2 (recalling that 
in (|5.5p we chose i = 1) we can obtain the following inequality : 



(7 £ )' k w l ° 2 ' £ dx + J^)' k < v 1 '* dx < rf(C + l)w\{x) +c(l+^j, 



(5.14) 



where we also used the fact that HifflU 00 < wf (a;). Summing inequalities (|5 . 1 3[) and (|5.14p and 
choosing 77 > small enough the conclusion follows. □ 

Next lemma gives a control of the Hessians D 2 u\ and D 2 u 2 in the support of a 1 ' 5 and er 2,£ 
respectively. 
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Lemma 5.7. There exists a positive constant C, independent of e, such that 



sup / e\D 2 u e .\ 2 a^ e dx < C. 

3=1,1 JU 



Proof. The bound of the Hessian of D 2 u\ comes from identity (|5.12l) . together with Lemma [ 
The other bound can be obtained in an similar way. □ 

We can finally prove the analogous of Lemma 12.21 
Lemma 5.8. There exists a positive constant, independent of e, such that 

max < C, max 7 e (0f (x) ) < C. 

xeu xeu 

Proof. It will be enough to prove the second inequality, since the first one will follow by the 
definition of 7 e . If the maximum is attained at the boundary, then 

max 7 e (0f (x)) = max7 E (-4(i)) = 0. 

j=l,2 J .7=1,2 
xeU xedU 

Otherwise, let us assume that there exists x\ £ U such that 

maxmax 7 £ (^) = 7 e (6lf)(2:i) > 0, 7 £ (^(iei)) = 0. 

Since 7 e is increasing and 7 e (z) > if and only if z > 0, we also have max^, g jj(6^ (x)) — 0\ (x±) > 0. 
Evaluating the two equations in (|5.3p at x\ and subtracting the second one from the first one 

6l( Xl ) + 7 £ (^i(^i)) = eAwf ( Xl ) - eAu|(a;i) - #i(zi, Du\( Xl )) + H 2 (x u Du e 2 {x x )) - ifrfa) 

< sAipxixt) - H^x^DuKx!)) + H 2 {x 1 ,Du 2 (x 1 )) - ifa(xi) 

< ||A^i(-)||£- + ||ffi(-,I>uf(-))IU- + \\H 2 {-,Du e 2 (-))\\ L oo + ||ViO)IU- < C, 

where we the last inequality follows from Lemma 15.61 □ 
We now set for every e £ (0, 1) 

du £ n _ 
ul £ (x) := -^fO), x£U,j = 1,2. 

The next lemma gives a uniform bound for u\ E and u 2 £ , thus concluding the proof of Theorem l5.ll 

Lemma 5.9. There exists a positive constant C > such that 

C 

max 14, (x) | < 
xeu 
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Proof. If the above maximum is attained at the boundary, then 

max Ittf = max lid Jx)\ = 0, 

3=1,2 J ' 3=1)2 
x£U xEdU 

since u\ E = u\ e = on dU. Otherwise, assume that there exists x € U such that 

P^KeWI = K,e(X)|- 
3=1^2 J ' 

Differentiating (I5.3P w.r.t. e we have 

(1 + (Y)' \ 0f )ul e + D p Hi ■ Dul e - (YY \ 6 e u|, E + 7 | | fl .= eAu\ t£ + Au\ in CT, 



(5.15) 



(1 + (YY h)u% e + D p H 2 ■ Dul £ - (YY |0f < e + 7s e |e S = eA4 e + A«| in U. 

Let cr 1 ' 5 and <r 2,£ be the solutions to system (I5.4[) with i = 1 and £0 = Multiplying (I5.15l) i 
and (|5.15p ? by a 1,£ and (7 2,<r respectively, integrating by parts and adding up the two relations 
obtained we have 

2 / r r 
^ ' ' ' %\e' V 



Thus, 



4 )£ (*) = E(/^A^' £ <fe- f u 



cl,M\ < E (/ y I Am 3 £ I ctJ,£ dx + J v l% £ l*j I ^ £ dx 



At this point, the proof can be easily concluded by repeating what was done in Section [5] showing 
relations d2H2) - (EUS]) . □ 
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